Screening in gated bilayer graphene via variational calculus 
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We analyze the response of bilayer graphene to an external transverse electric field using a vari- 
ational method. A previous attempt to do so in a recent paper by Falkovsky [Phys. Rev. B 80, 
1134f3 (2009)] is shown to be flawed. Our calculation reaffirms the original results obtained by one 
of us [E. McCann, Phys. Rev. B 74, 161403(R) (2006)] by a different method. Finally, we generalize 
these original results to describe a dual-gated bilayer graphene device. 
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The physics of monolayer and bilayer graphene has 
been a subject of much interest recently—' A unique fea- 
ture of bilayer graphene (BLG) is its tunable band struc- 
ture: its bandgap depends on the external electric field, 
which can be controlled by doping or gating. This effect 
was first analyzed theoreticallySSl and recently studied 
experimentally."^ ^ 

The subject of this paper is the value of the gap at 
the Brillouin zone corners, which we denote by 2|C/|. As 
pointed out previously,^ ^ 2U coincides with the electro- 
static energy difference per electron in the two layers. 
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where is the component of the electric field inside 
the BLG directed from the bottom to the top layer and 
dm is the interlayer spacing. This field depends not only 
on the external gate voltages but also on the induced 
electron densities rit and nt of the layers, see Fig. [T] For 
a given chemical potential fi, these densities, and thus the 
total electron density of BLG n = nt + rib are nonlinear 
functions of U. I n general, they can be calculated only 
numerically.l^'^E^'^Sl However, in a range \U\ <^ ^ 71 
an asymptotic analytical formula for the interlayer bias 
was derived by one of u^ 



N -2N, 



dt 



A 



Mc 



\N\-nn\N\ 



N : 



(2) 



The derivation was done within the commonly used ap- 
proximation that neglects certain small electronic struc- 
ture parameters 73, 74, and A' (for their physical mean- 
ing and a discussion of their numerical values, see 
Ref. ,20,) The parameters that are retained include the 
interlayer hopping energy 71 = 0.41 eV and the Fermi ve- 
locity of monolayer graphene v — 1.0 x 10^ cm/s, which 
define the characteristic density scale in the problem 
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The remaining notations in Eq. ^ are 
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FIG. 1. (Color online) Device schematic. The BLG is shown 
as two horizontal lines in the middle. The thicker line at the 
bottom represents the gate. The arrows show the electric field 
in the system and the crosses depict positive charge density 
Udt on (or on and above) the top layer of the BLG. 



ndt on the top layer (or above it, see below) and 
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which is the dimensionless strength of the interlayer 
screening, with Km ~ 1 being the effective dielectric con- 
stant of the medium between the layers. The screening 
is particularly significant in the narrow gap regime. If 
Ndt = 0, it is realized in the limit n ^ Q. According 
to Eq. ([2| , the screened field Em is suppressed compared 
to the external field Ei, by a divergent logarithmic factor 
[see Fig. [1] and also Eq. (37) below]. As a result, 2U has 
a superlinear dependence on n. 

Formula ^ was recently challenged by Falkovsky.l^ 
Under the same conditions and assumptions, he obtained 
the simple linear law: 



2U = C7iA^ , 
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The goal of this paper is to show that Ref. contains 
egregious mistakes. Once they are corrected, Eq. ^ is 
recovered. 

Some of these mistakes are apparent on physical 
grounds. One key formula, Eq. (8) of Ref. [21] states 
(in our notations) that the BLG polarization 



which is the scaled background density of positive charge 



An = nt — nb 
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is positive for U, ji > 0, i.e., the BLG is polarized against, 
not along the electric field. In principle, such a phe- 
nomenon can arise if a system has a negative compress- 
ibility'^^''^ due to exchange and correlations (see the dis- 
cussion at the end). However, neither exchange nor cor- 
relation effects are included in the theory of Ref. [2TJ 

How the author of Ref. f5T is able to eventually arrive at 
Eq. ^ is also unclear to us. Even if we take for granted 
the basic Eqs. (8)-(10) of Ref. and dutifully follow 
the suggested steps of the derivation, we are still unable 
to reproduce the later equations. For example, on the 
left-hand side of Eq. (13) of Ref. we get the top layer 
density instead of the bottom one stated there, and on 
the right-hand side we get an extra minus sign. With 
none of the choices for the sign or the layer can we get 
from Eq. (13) to Eq. (14) of Ref . [HI 

To settle the matter decisively we carry out below our 
own derivation using the method proposed by Falkovsky 
in Ref. [211 Its idea is to treat U and /i as variational 
parameters and get their actual values by minimizing the 
total energy density E of the system at a fixed n. Note 
that fj, and U are admissible variational parameters if the 
Jacobian of the transformation from rit and ni, to fi and 
U is non-degenerate: 
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This condition can be verified from the final results. We 
show that the variational method leads back to Eq. ^ 
and that all the intermediate results make physical sense. 

Let us proceed. In the mean-field approximation E is 
the sum of the kinetic and the Hartree interaction terms: 



E — -Ekin + E^ , 

The Hartree term is straightforward: 



(9) 



En = ]^e^Cb{n - n^J 



[nt-ndtf, (10) 



which can be understood as the energy of two parallel- 
plate capacitors in series. Here — ATTHb/db is the ca- 
pacitance per unit area between the BLG and the bot- 
tom gate, Kb is the dielectric constant below the BLG, 
and n^b is the density of additional positive background 
charge charge on the bottom layer (not shown in Fig. [ij . 



For = 0, Eq. pl)|) agrees with Eq. (10) in Ref. 



The kinetic term i?kin requires a little care. The author 
of Ref. [2n seems to assume that it coincides with the sum 
of the occupied single-particle energies 



(11) 



where a is a short-hand notation for all quantum num- 
bers, Q{x) is the unit step-function, and A is the area of 
the system. In fact, each energy Ca, which is an eigen- 
value of the BLG Hamiltonian H = 7Jkin+$ also contains 
a potential term: 
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Here $ stands for the electrostatic energy that is equal 
to zLU in the top (bottom) layer. Lumping together Eg 
and En is incorrect as it leads to double-counting of the 
interaction energy. Instead, the proper formula is 

E^in = ;^(a|i7ki„|«) - - UAn . (13) 
Next, using the Hellmann-Feynman theorem, we get 



An = > ( a — 

Note also the useful relations 

n = nt + Ub = 1 
Therefore, 
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The correct variational functional to minimize is 

a 
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where ^ is a Lagrange multiplier. The sought minimum 
is determined by the equations dft/dfi = dil/dU — 0. 
Using Eqs. ( 16 ) and ( 15 ), one can bring them to the form 
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In view of Eq. ([8 1 , the system of linear equations for C 
and D has only the trivial solution C — D — and so 



2U = ■i^Te^{d„Jn„^){nt - n^t) ■ 



(21) 



Of course, this equation can be written at the outset be- 
cause it follows from Gauss's law. Indeed, it is the start- 
ing equation [Eq. (2)] of Ref. ffl The closest to Gauss's 
law in Ref. 21 is Eq. (13) of that paper but it has the 
problems described above. 

Let us now compute n^ by the present method to 
demonstrate where Ref. contains another mistake. 
The total polarization An is the sum over all four bands 
of BLG 
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We use the signature ai, (72 to label the bands as follows: 
CTi = ± distinguishes the conduction and valence bands 
and (72 — ± labels the outer (inner) bands. The band 
dispersions are given b}0 



([/, fc) = ai W ^ + £2 (fc) + [/2 + a2i?e 



^+eHk){jl+4U^)^ 
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where k is the deviation of the quasimomentum from the 
nearest Brihouin zone corner and e(fe) = hvk. Each band 
has 3 = 4 fold degeneracy due to spin and vaUey. 

For simpHcity, we consider the case where the Fermi 
surface is singly-connected and includes only the states 
of the ai,<T2 — band, so that n = g{TTkp)/{2Tr)^ = 
kp/n where kp is the Fermi momentum and 



+ £!i|£i ^ + ^2^2 . (25) 



As explained above, if U is positive, then An should be 
nega tive. In Ref. [21] however only the conduction band 
term,lil2il 



An+,_ = In - 

71 \ U 




(26) 



is included while the much larger negative contribution 
An_ = An_^+ + A?t._._ of the completely filled valence 
bands is neglected. This important term is given by 
[cf. Eq. Jm])] 



d 
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The above integral can be done exactly but we need only 
its leading asymptotic form: 
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which is in agreement with Ref. |31 Combined with 
Eqs. (251 and (26), it yields the result of the correct sign. 
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Substituting it into nt = (n + An)/2 and using Eqs. (|4]), 
([5]), and (21), we obtain the equation for the function 
U{N): 



2[/(7V)~7i- 



N - 2Ndt 



(30) 



A 



Mc 



In 



This formula is actually valid for arbitrary signs of U and 
N . In the region 



2C/(0)/7i < TV < 1 



(31) 



it coincides with Eq. ^ within the accuracy of this cal- 
culation. Here 
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is the gap at iV = 0. Note that the gap vanishes only at 
one density, N = 2Ndt, i.e., n = 2nrft, at which 

nt = Ub = Udt ■ (33) 
The gap remains nonzero at all other n, in disagreement 
with the claim made in Ref. 21 but in agreement with 
the numerical results of Ref. ,4, 

As mentioned in the beginning of this paper, the back- 
ground positive charge density n^t does not have to re- 
side directly on the top layer. Equation (30) remains un- 



changed if some of this charge is located above the BLG, 
as sketched in Fig.[T] For example, for a dual-gated BLG 

^{eVt-fi + U), 
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where the first term is the fixed donor density on the top 
layer and the second term is the (tunable) charge density 
on the top gate. Parameter ct is the capacitance between 
the top gate and the BLG per unit area. Similarly, the 
total density is given by 



C* - Cf, 
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where Vb is the voltage on the bottom gate. The chemical 
potential /i enters these equations because the measured 
"voltages" Vb and Vt of the two gates are not simply 
the electrostatic but instead the electrochemical potential 
differences between these gates and the BLG, cf. Refs.l^ 
andL24_, However for distant gates, 

one can in the first approximation neglect U and /i com- 
pared to eVb and eVt , leading to the simplified equations 

eE^dm 
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(38) 
(39) 



which may be useful in experimental practice. Note how- 
ever that unavoidable disorder creates additional correc- 
tions to these expressions, l^fil which can be as large as 
30%. 

In closing, we comment on electron-electron exchange 
and correlations. These man y-b o dy e ffect s have been pre- 
dicted to generate a band gajji^I^llEZl and the correspond- 
ing polarization An 7^ even in the absence of the exter- 
nal electric field. Although neglected in our mean-field 
theory, exchange and correlations can be incorporated by 
introducing a suitable self-energy difference between the 
layers, which is a function of N and U. In fact, we already 
have a similar parameter in our formalism — this is our 
Ndt ■ We intend to investigate this interesting problem in 
a future work. 
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